On the Incommensurate Phase of Pure and Doped Spin-Peierls System CuGeC>3 
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Phases and phase transitions in pure and doped spin-Peierls system CuGeOa are considered on 
the basis of a Landau-theory. In particular we discuss the critical behaviour, the soliton width 
and the low temperature specific heat of the incommensurate phase. We show, that dilution leads 
always to the destruction of long range order in this phase, which is replaced by an algebraic decay 
of correlations if the disorder is weak. 

PACS numbers: 64.70.Rh, 75.10.Nr, 75.30Fv, 75.40.-s 
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The spin-Peierls (SP) transition is the classic instabil- 
ity of one dimensional quantum spin-half antiferromag- 
netic chains due to the coupling of the spins with the 
lattice. A rigid Hcisenberg chain has a nonmagnetic uni- 
form ground state with a gapless fermionic excitation 
spectrum jjj. This can be seen most easily by using 
the Jordan- Wigner transformation, which maps the spins 
onto (strongly) interacting pseudo-fermions |Q. Due to 
the coupling to the lattice the system can lower its energy 
by the standard Peierls mechanism: lattice distortions 
freeze in at a wave vector 2kp which leads simultane- 
ously to the opening of a gap at the Fermi-level in the 
fermionic spectrum such that the energies of all occu- 
pied fermionic states decrease [|| . In zero magnetic field 
the free-fermion band is half-filled with Fermi-wave vec- 
tor Uf = sj, which corresponds to a dimerization of the 
chain. A non-zero magnetic field lowers the Fermi-level 
||, but Umklapp processes still favor the distortion at 
ir/a until a critical field strength Hj is reached, at which 
a transition to an incommensurate (I) phase with mod- 
ulation vector \2kp — q s \ sets in. In the I-phase a new 
(empty) band appears in the middle of the gap of the 
fermionic spectrum. Thus, spin excitations still exhibit a 
gap which is however smaller than the gap of the dimer- 
ized (D) phase. The above picture follows from theories 
obtained for free or weakly interacting pseudo fermions, 
in which phonon dynamics were essentially ignored ^],[| . 
There, the SP transition is the result of the freezing of a 
(classical) phonon mode due to further downwards renor- 
malization of the phonon frequency by the spin-phonon 
interaction ^[(^]. This scenario is supported by the ex- 
perimental data of organic SP systems 0. However, it 
does not seem to apply in all respects to the transition 
found recently in the inorganic SP substance CuGeC>3 [|| . 
Though not devoid of controversy, there is now a wealth 
of well-accepted results for CuGeOs, which shows two SP 
transitions |] [l3|,|l5|,|l6| . The SP transition from a disor- 
dered, uniform (U) to a D phase at 14.3 K in zero field 
is shifted slightly to lower T if the field increases until a 
Lifshitz point at T w 11.3 K and H w 12.5 T is reached, 
where the transition to an I phase sets in. Some exper- 



imental results which are not explained by the existing 
theories, are: (i) no soft phonon has been observed so far 
§,|, (ii) a (Peierls) gap in the D phase is observed in low 
temperature specific heat measurements, but not in the 
I phase, for which a Debye-like T 3 -law has been found 
with an amplitude much larger than the background (lat- 
tice) contribution [ pof , (iii) solitons, which are supposed 
to produce the modulation in the I phase are broad |Tl| ] 
in comparison to the sharp Sine-Gordon like solitons pre- 
dicted by mean- field like calculations S, (iv) already a 
small amount of doping leads to a strong reduction of the 
SP temperature T$p |l6| and a drastic suppression of the 
anomalies at the UI transition 

Since the phonon energies are always large compared 
with the magnetic ones, the applicability of the adiabatic 
approximations has been questioned pl| . Khomskii et 
al. |l^] developed a soliton picture of the SP transition 
in CuGeOa, which resembles somewhat structural order- 
disorder transitions p3]. No soft mode phonon is ex- 
pected, but the SP transition corresponds to deconfinc- 
ment of solitons, which are bound to pairs below Tsp. 
These solitons are simultaneously magnetic and struc- 
tural excitations: they carry spin 1/2 and are domain 
walls between the two groundstates of the dimerized lat- 
tice. 

It is the aim of the present paper to explain the prop- 
erties (ii)-(iv) by a pure phenomenological approach, 
which avoids delicate approximations in the coupled spin 
phonon system: The T 3 -law of the specific heat in the I 
phase is explained quantitatively by phason fluctuations. 
It is argued that broad solitons are fingerprints of the 
type II lock-in transition which occurs in SP systems like 
CuGeOa. Finally, we show that dilution leads to com- 
plete destruction of long range order in the I phase. 

Incommensurate phases are classified according to the 
existence of an inversion symmetry for the structural 
transition in question |21|. In case there is an inver- 
sion symmetry for the Hamiltonian, as for CuGeOs, 
first derivatives of the order parameter (Lifshitz invari- 
ants) do not exist. Indeed, for CuGeOs the uniform 
high temperature orthorhombic structure, space group 
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Pbmm, changes below Tsp to a dimerized structure, 
space group Bbcm, with distortion wave vector (55O) 
(established from X-ray and neutron diffraction experi- 
ments ||). Standard group theoretic arguments based 
on the symmetries and the invariant group of the dis- 
tortion vector in the Brillouin zone |22| show, that the 
transition is described by four non-equivalent, one dimen- 
sional irreducible representations |2^]. It is very likely, 
that only one of these four representations corresponds 
to the primary order parameter, which is real and can 
be considered to be proportional to the displacement of 
the copper ions. The other three may occur as secondary 
order parameters. In fact, neutron scattering data indi- 
cates that at least two normal modes are necessary to 
explain the displacement pattern of the D phase fl. A 
simple transformation of reversing the displacements on 
one sublattice helps us in getting an ordered state with 
zero wavevector for the D phase. We take this trans- 
formed and coarse grained displacements to be propor- 
tional to the order parameter field i/j(x). In general, ^( x ) 
will also include contributions from the magnetic degrees 
of freedom. 

The Landau Hamiltonian is that of an anisotropic Ising 
model HII = / d d x/i{V>(x)} , where the Hamil- 

tonian density is (for i = x,y, z) 

+ ^») 2 + |^ + |v 2 (^) 2 (i) 

Terms involving d z ip d z ip, though not shown explicitly, 
may also appear We have included higher order 

terms in order to stabilize the system for the case when 
one of the Cj coefficients (here c z ) or u becomes nega- 
tive for sufficiently strong magnetic field. A negative Cj 
signals the transition to the I phase. The parameters r 
and Ci are taken as analytic functions of T and H with 
r = r {T- T SP (H)), and c z (T,H) = c (JJi(T) - H). 
Microscopic treatments || and experiments suggest [[lOj , 
that u also decreases considerably with increasing mag- 
netic field. 

A mean field analysis (MFA) of the phase diagram re- 
quires treating the Hamiltonian Eq. (^) as a free energy 
minimized with respect to ip. Let us first assume, that 
u and d remain positive everywhere in the H T plane. 
Then, we can ignore the last two terms in ([!]) . If Cj > for 
all i, the minimum of the free energy occurs for wave vec- 
tor k = 0, while a nonzero k vector is possible if c z < 0. 
The mean field phase boundaries are given for UD: r = 0, 
c z > 0, UI: r = c 2 z /2d and DI: r = \{y/3/2 - \)- x (?Jd 
p4j . The spontaneous wave vector in the I phase is given 
by q 2 = —c z /d. The DI transition is first order while 
the other two are second order, in agreement with exper- 
iments. 

Fluctuations do make subtle changes in the phase di- 
agram but the overall features remain the same. The 



critical behavior of the UD transition is Ising-like but 
the UI transition is XY-like. This difference originates 
from the fact, that in the I phase the order parameter 
condenses at k = ±q s and has consequently two compo- 
nents A(x) = (Ai(x), Aa(x)). These are related to t/>(x) 
by tp(x) = \/2 [Ai(x) cos(q s x) + A 2 (x) sin(q s x)]. In the 
I phase, the Landau functional thus can be written as 

h{A} = |(r - c 2 z /2d)A 2 + 2\c z \(d z A) 2 

+ c x (d x A) 2 +c y (d y A) 2 +lu(A 2 ) 2 (2) 

Since the number of degrees of freedom of the system 
cannot change when going from the D to the I phase, it 
is clear that Eq. (||) is valid only for fluctuations of A(x) 
with wavelength long compared to q~ , i.e. as long as we 
are away from the Lifshitz point. 

If one approaches the ordered phase along the line 
c z (T,H) — 0, one observes so called Lifshitz critical be- 
haviour, which follows from a change of the dispersion 
relation to A-^ = r + c x k 2 + c y k 2 + dk z . Note that at 
the Lifshitz critical point the conventional hyperscaling is 
changed to + (d— l)v±_ = 2 — a where v\\ = v^/2 = 0.31 
are the correlation length exponents parallel and perpen- 
dicular to the ^-direction p4| . Approaching the D or I 
phase, respectively, from the U phase on a line parallel 
to that given by c z (T, H) = 0, at first Lifshitz type criti- 
cal behaviour will be observed before the region of Ising - 
or XY-type critical, respectively, behaviour is asymptot- 
ically reached. 

Considering the DI transition, fluctuation effects are 
expected to be less important, because it is first order 
in MFA. A refined MFA has been worked out by Bruce, 
Cowley and Murray |2l| for this case, who found that 
in the I phase the order parameter can be described by 
a multiplane-wave Ansatz ^(x) = S °™ cos (mqzz) with 
m = 1, 3, 5..., which is rapidly converging. For example, 
the ratio |a 3 /ai| ~ 0.035 close to the transition In 
this sense the system shows broad domain walls. Also in 
this refined theory the transition remains first order. 

Above we assumed u to be positive even for large field 
values. In the opposite case, the transition to the U phase 
might become first order. Some mean-field theories 
predict very special relations between the coefficients of 
the Landau-expansion, i.e. 

u/c z = const. , w — 3du 2 /4c 2 . , e = 5du/2c z . (3) 

If these are fulfilled, the DI transition may become con- 
tinuous, at least close to the Lifshitz point j25|. In- 
deed, for this very particular relation of the coefficients of 
the Landau-expansion (^) , the ground state solutions are 
the Jacobian elliptic functions V's(z) = sin((/) s (z)/2) = 
ipo sn(z/fc£ s , k), where £ s is a bare correlation length (ex- 
pressed by c z , d and r) and k is the modulus of this 
function ||^6]]. Note, that 4> s {z), which is related to the 
spin-density, obeys the Sine-Gordon equation. In this 
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case solitons are sharp in the sense that the separation of 
domain walls diverges by approaching the D phase. How- 
ever, from a symmetry point of view, which we adopt 
here, we do not see a deeper reason, why the relations 
(|J) should be fulfilled in general by an exact microscopic 
theory. In fact, these relations were obtained using the 
adiabatic approximation. Consequently one has to ex- 
pect, that in general the w- and e-terms in (|l|) do exist, 
but violate the relations (||). These terms will change the 
modulation amplitude ratio |a3/ai| to larger values, but 
without reaching the sharp soliton limit. Thus the DI 
transition is expected to remain first order, as found also 
experimentally for CuGeC>3 

Although the validity of Landau-theory is essentially 
restricted to the region close to the transition, one should 
expect that it can be used to understand at least qual- 
itatively the low temperature specific heat data. For 
this purpose, we have to determine the low-lying exci- 
tations of the ordered structure. These can be found by 
adding the kinetic energy term J d d x|-0 2 (x) to the GL- 
Hamiltonian, where the mass density p will have contri- 
butions both from the magnetic and the lattice degrees 
of freedom. We will further assume, that ip obeys Bose- 
statistics. Using the saddle point approximation to de- 
termine the equilibrium value of ip 1 one obtains in the 
D phase w 2 (k) = j(2\r\ + Cjfcf) for the frequency of the 
harmonic excitations of the order parameter field. In 
the D phase, were the order parameter is real, we iden- 
tify Eg — h(2\r\/ p) 1 / 2 with the gap which is found in 
the low-T specific heat. In the I phase in addition to the 
massive amplitude mode a gapless phason mode with fre- 
quency w 2 (k) = i {c x kl + Cykl + 2\c z \k 2 z ) appears 0, 
which will dominate the specific heat 
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Here we have introduced £ = (Coz£oy£oz) where 
£oi = c i/ r oTo and used T « T /2 to express p by E g . 
Thus, the phason mode delivers a T 3 -contribution to the 
low-temperature specific heat, in addition to that from 
acoustic phonons 

Next we extend our analysis to the quenched disor- 
dered case, e. g. random substitutions of Cu by Zn or Ni 
and/or Ge by Si in CuGe03. Such substitutions change 
the various interactions locally but do not break the sym- 
metry of the displacements in favor of a particular dimer- 
ization. Therefore, the effects of these random substitu- 
tions can be modeled by randomness in the coefficients of 
the original Landau Hamiltonian without any symmetry 
breaking term. Little reflection shows, that the main ef- 
fect will come from a randomness <5r(x) in r p8[. In the 
D phase, this leads to a decrease of Tsp, as was shown 
microscopically by Khomskii et al. [M. Moreover, the 
critical behaviour will be changed to that of the diluted 
Ising model p9[. 



The effect of disorder is even more severe in the I phase. 
This can be seen easily by rewriting ^(x) as ip(x.) — 
v / 2A(x) cos (q s x + 0(x)). With 5r(x) = k£ ( <*(x - x 4 ) 
the disorder term can now be written as 



|^A 2 (x,) cos [2 (0(xi) +q s x 4 



(5) 



The random impurity positions xj lead to a random 
phase cti = a(x,-) = 2q s Xi (mod 2n) which is equally dis- 
tributed between and 2n. It is well known that such a 
random anisotropy term destroys the translational long 
range order of the I phase p0| |. However, the phase- 
phase correlation function diverges only logarithmically 
§l| ((0(x) -0(O)) 2 ) = hx{x/L L ). Here the overbar 
denotes the disorder average. The Larkin-lcngth is 
related to the strength of the disorder, a rough estimate 
is 



L L » 2tt 3 [$/{dlnT SP /dn imp )] 
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where nj mp denotes the concentration of the impurities. 
Because of the logarithmic divergence of the phase fluc- 
tuations, there is however quasi-long range order of the 
order parameter correlation function 



A(x)A(0)) « W£i(*i/£o,i) 2 /££ 



-tt 2 /36 
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Despite of the loss of true long range order, the system 
will however still show Bragg peaks of finite width, as 
follows from the Fourier transform of (j^). In deriving 
these results we have neglected vortex-ring excitations. 
It has been argued recently, that these can indeed be 
neglected for sufficiently weak disorder and low temper- 
atures |32|| . At elevated temperatures or larger dilution 
their condensation triggers the transition to the disor- 
dered phase. The type of this transition is presently un- 
known. 

We briefly apply the results obtained so far to CuGe03. 

Fixing the T = value of the order parameter at 

ipo = 1, we have r T a = u . From the mean-field jump 

of the specific heat ACmfa = r oT Q /2u = u Q /2T w 

22.7 ^^3 in zero magnetic field JlO|| , we get u — 

650-2^ which gives the correct size of the critical re- 
cm" 3 ° 

gion. Since ACmfa decreases for increasing field and 
is reduced approximately by a factor 4.6 when reach- 
ing the Lifshitz point, u is reduced correspondingly to 
about ul — 112^3, but still positive. Defining the 
Ginzburg critical region tq = \Tq — Tq \ /Tq as the re- 
gion, in which the first fluctuation correction to the spe- 
cific heat becomes larger than ACmfa, this gives for 
zero field, ^q jX = 0.12 nm, £o,y = 0.36 nm and £0,2 = 
0.69 nm p4[ , with geometric mean £ = 0.31 nm, tqj ~ 
(fc_B7o/87TM^Q) 2 w 0.16 - larger ACmfa diminishes tq 
correspondingly. For the XY-transition far from the Lif- 
shitz point we get tq.xy ~ 0.32 at a magnetic field where 
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£,z{H,T = T SP ) w f ,*- At the Lifshitz point the critical 
region is given by t g , l w t^ 3 (£ ,2U /\/2C"l) 4/3 « 0.7 
where Co = (d/2roTo) 1 ' 4 ~ 1.2 nm. The critical exponent 
(3 changes from /3j ~ 0.325 for H < H L to (3 L ~ 0.15. .0.18 
for H ^ H L and then to /5 XY ^ 0.346 for H > H L , in 
agreement with the experimental observation |D| . 

From the low-temperature specific heat in the D phase 
one finds £ s 23K which gives with Eq. (g) 

for the phason specific heat /3 p hason ~ 1.3 K ^ ol in the I 
phase, which has to be compared with the experimental 
value of 1.4 K "^ ol PJ- This good agreement is possibly 
to some degree accidental, since the magnetic field de- 
pendences of the various parameters have not been taken 
into account carefully. But at least the order of magni- 
tude should be right. 



For the Larkin length we obtain with x = n 



/2 



imp u uc 

for the concentration of the Zn-atoms (v uc denotes the 
volume per unit cell) and assuming a linear dependence 
of T SP {x) on x with d\nT SP (x)/dx « 14 |§, for x » 
0.04 : L L rj 1.2 nm and for x « 0.07 : L L « 0.7nm. The 
data of Kiryukhin et al. [p"i"| was fitted with an exponen- 
tial decay of correlations with an anisotropic correlation 
length £ of order 10 nm. It would be interesting to check, 
whether their data can also be fitted by a power law (]?]). 
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